Cutting-Decimation Renormalization for diffusive and vibrational dynamics on fractals 
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Recently, we pointed out that on a class on non exactly decimable fractals two different parameters 
are required to describe diffusive and vibrational dynamics. This phenomenon we call dynamical 
dimension splitting is related to the lack of exact decimation invariance for these structures, which 
turn out to be invariant under a more complex cutting-decimation transform. In this paper we study 
in details the dynamical dimension splitting on these fractals analyzing the mathematical properties 
of the cutting-decimation transform. Our results clarify how the splitting arises from the cutting 
transform and show that the dynamical dimension degeneration is a very peculiar consequence of 
exact declinability. 
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I. INTRODUCTION 



Fractal dynamics is a field of primary importance in the study of physical phenomena in real non-crystalline systems. 
The analytical study of dynamical properties of fractals requires the introduction of simplified model structures, which 
are supposed to have the same universal properties of the real ones. However, as a matter of fact, analytical results are 
usually restricted to exactly decimable fractals, where one can apply the powerful techniques of exact renormalization 
group. Exactly decimable fractals are indeed very peculiar structures, characterized by strong restrictions on their 
topology, which are far from being general and representative of all fractals. Nevertheless, due to the availability of 
analytical results, their behavior was usually considered as typical. This is the case for the spectral dimension of a 
fractal. Such a parameter, the spectral dimension d, was defined by the asymptotic law 0]: 

pM ~ J- 1 (l) 

where p(uS) is the density of harmonic vibrational modes with frequency u>, or by: 

P«(t)~H (2) 

where Pu(t) is the probability of returning to the starting site i after t steps ( t — > oo ) for a random walker and 
the exponent d/2 is independent of the starting point In the original definition, ([j]) and were supposed to be 
equivalent by scaling arguments and in all calculations made on exactly decimable fractals they have been found to be 
always equal. However it has been shown jl that from a mathematical point of view it is not possible to conclude that 
the two asymptotic behaviors of (0) and (|2|) are equivalent. Therefore, one should distinguish between the diffusive 
spectral dimension defined by (g) and the vibrational spectral dimension defined by ([!]), the former being relevant for 
local quantities and the latter for bulk or average quantities. In the following we will show that the coincidence of 
the two spectral dimensions is typical of exactly decimable fractals while more general structures present dynamical 
dimension splitting Q]. Indeed for exactly decimable fractals p{u>) and Pu{t) have the same transformation properties 
under renormalization and this leads to the so called spectral dimension degeneration. Here we will consider the 
case of a class of non exactly decimable fractals (called NTd) showing that (|l|) and (§) can be calculated applying a 
generalized renormalization transformation, different from the usual decimation and leaving the structure invariant. 
The key point is that here p{u>) and Pu(t) transform according to different laws under this new renormalization 
procedure. This leads to independent asymptotic behavior for these quantities and requires the definition of two 
distinct spectral dimensions we will call do and dy, D staying for diffusive and V for vibrational. 



II. HARMONIC OSCILLATIONS AND RANDOM WALKS ON GENERIC NETWORKS 

The harmonic oscillations of a generic network of masses m connected by springs of elastic constant K are described 
by the equations of motion for the displacements x% of each mass from its equilibrium position: 

d 2 

where the sum runs over the nearest neighbors of point i. 

Equations (W can be Fourier transformed with respect to the time giving: 



where u> 2 = K/m. 

The equations describing random walks and harmonic oscillations on a network are formally similar. Let us consider 
a discrete time random walk (the so called blind ant problem) on a network; the master equation for the probability 
of being at site i after t steps for a random walker starting from an origin site O at time 0, is: 

Poiit + 1) - P 0i (t) = l/z max J2( p Ooit) - Poi{t)) 

P Oi (0) = S i0 (5) 
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is the maximum coordination number of the network (i.e. the maximum number of nearest neighbors of 
a point), that is assumed to be finite. Applying the discrete Laplace transform with respect to the time defined by 



£i(e) = £(l + e)-%-(i) (6) 

and taking e — > 0, the system (|J) can be written in a form similar to (Q): 

z max eP i(e) = ^2(?Oj(e) - Poi{e)) + 5 0l - (7) 



After the substitutions —z max e — > uj 2 /uj 2 . — 7 and Poiif) — > %i, equations (|j) and (Q) look very similar. However they 
present a fundamental difference consisting in the term Sen, arising from the initial condition in (^). In addition, the 
system (^) is homogeneous and defines an eigenvalue problem which has infinite solutions (all the normal modes of 
the graph). The density of eigenvalues for lu — > depends on dy and it is not a direct solution of the system, but has 
to be separately calculated after solving the whole system. On the other hand, the system (f7|) is inhomogeneous and 
corresponds to a Cauchy problem, which has only one solution. The behavior of such a solution for e — > depends on 
do and, in principle, has no direct relation with the density of vibrational modes and with dy. A relation between 
the random walk probabilities Pa (t) and dy does indeed exist but it involves the average of the Pa (t) over all points 
of the graph [||: 

1 N 

P(t) = lim — V Pii(t) ~ r dv ' 2 for t^oo (8) 

N->oo N ^ 
i=l 

Relation (|J) does not imply that do = dy since the average P(t) has in general a different asymptotic behavior from 
each Pu(t). 



III. RANDOM WALKS AND HARMONIC OSCILLATIONS ON EXACTLY DECIMABLE FRACTALS 

Exactly decimable fractals are a restricted class of self similar structures (i.e. not all self similar structures are 
exactly decimable) which are geometrically invariant under site decimation. Examples of exactly decimable fractals 
are the Sierpinski Gasket 0, the T— fractal the branched Koch curves Q and so on (Fig.l). The solution of 
both the random walks and the harmonic oscillations problems can be obtained by standard renormalization group 
calculations based on a real space decimation procedure S. 

A geometrical structure is decimation invariant if it is possible to eliminate a subset of points (and all the bonds 
connecting these points) obtaining a network with the same geometry of the starting one. From a mathematical point 
of view this corresponds to the possibility of eliminating by substitution a set of equations from system (||) or (0) 
obtaining a system which is similar to the initial one after a suitable redefinition of the coupling parameter 7. If we 
consider for example a T— fractal the decimation procedure consists in transforming each "T" made of three bonds 
in a simple bond connecting two points. As it can be easily verified this operation does not change the geometry 
of the network but requires a redefinition of the coupling parameter 7. In general, after a decimation step, 7 splits 
into a finite number of different couplings 7^, ji = 1, n on geometrically inequivalent points. In the case of the 
T— fractal there are two kinds of inequivalent points: points having one nearest neighbor and points with three nearest 
neighbors. This suggests to distinguish between a coupling 71 and a coupling 73 to be used in equations (ji|) or (Q) 
where point i has respectively one or three nearest neighbors. Before the decimation we have obviously 71 = 73 = 7 
but this distinction is useful to put into evidence the forthcoming splitting. 

The splitting of 7 in an at most finite number of couplings is a necessary condition for exact decimability. If this 
condition is fulfilled, the decimation transform can be iterated and linearized near the fixed point 7^ = 0. After the 
linearization the transformation laws for the T— fractal are 

7i -* 7i = 3 7i + 173 
73 -> 73 = 3 7i + 573 (9) 
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The linearized decimation transform can be represented by a matrix D acting on a vector with components 7^, fi = 
1, n so that: 

^ = E^> ( 10 ) 

In our particular case: 

D = (l 5) 

Decomposing 7 M on the basis of eigenvectors of D, as the number of decimation steps goes to 00, 7^ tends to the 
eigenvector corresponding to the largest eigenvalue of D. Therefore the largest eigenvalue of D, a 2 , determines the 
transformation laws of coupling parameters. 

Following these steps in the case of random walks one finds for the parameter e the transformation: 

e^e'(e)~a 2 e (12) 

where e is now the projection of the vector e M on the largest eigenvalue direction. The presence of the term Sio in 
(0) requires a redefinition of the quantities Pij{c) to assure that, even after the decimation, the initial condition will 
correspond to the probability of being in a fixed site equal to 1. One introduces a new parameter c and writes the 
transformation law for Py(e) as: 

P^^P^e')^-^) (13) 
The diffusive spectral dimension d,D is obtained using the relation: 

Poo(e) ~ e 1 *' 2 - 1 (14) 

which holds only for do < 2. As we will discuss later, this is always the case for exactly decimable fractals. Now one 
can rewrite ( |l3|) as: 

P<i(e)~c%(0 (15) 

and observe that since our graphs are infinite Poo{t) and Poo( e ') re f er to the same structure and they must have 
the same functional form. This gives: 

e 7 D /2-l = c(a V)^/2-l ( 16 ) 

so that: 

~dn = 2^ (17) 
log or 

As for harmonic oscillations, the invariance of the network under the decimation procedure gives for u> 2 the same 
transformation law obtained for e. The analogous of @ is here: 

lo — > u> (u>) r~j aui (18) 

The relation between the density of vibrational modes p(o>) and the new p'(w') is given by: 

p{uj)duj — p'(uj')duj' (19) 

Since in the decimation procedure the initial set of equations (|J) is reduced of a factor r, condition ( |l9|) leads to the 
relation: 

p(w) -> pV) - -pM (20) 
a 



From (EFj) it follows: 
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p(w) = -p'(co') (21) 
r 

so that: 

u/v- 1 = -{aujfv- 1 (22) 
r 

and: 

dv = ^ (23) 
log a 

Comparing ([l?]) and (^3|) one realizes that dy = dp if the decimation ratio r is given by: 

r = a 2 /c (24) 

This can be shown to be the case for exactly decimable fractals, using results obtained for the Gaussian model on 
exactly decimable fractals. The Gaussian model on a graph is defined by the Hamiltonian: 

The autocorrelation functions of the model are related to the generating functions Pa(e) of random walks by: 

Pii(e) =< fcfc >{m 2 i } ~Y~^~e ^ 

where Zi is the coordination number of site i and the masses mf are given by mf = Zie/ (1 — e). Now it can be shown 
that a decimation procedure implies the following scaling relations for the coupling J and the masses mi: 

J — > aJ mi — > Prrii (27) 

where the parameters a and (3 are related to the expression of the diffusive spectral dimension by: 

T 21og5 , x 

log/3/Q 

Now the scaling relations ( p7j ) can be rewritten as scaling relations for the fields <pi and the masses mf. 

P 

(pi — > Va^i m i - * ~ m i (29) 



In terms of random walks the transformation on the masses in (|29|) does not affect the asymptotic behavior of P«(e) 
H while the first of (|29|), thanks to relation (|2^), implies that Pu(e) — * aPu(e) so that we can identify a and 1/c. 
Since the parameter /3 is the decimation ratio r, from (28), ( ]23| ) and (^) equality (|ij) follows. 

Notice also that since it is always f3 > a ^ 1 j|, exactly decimable fractals have do = dy < 2. From (|2~i|), the 
dimensional degeneration = dv can be obtained knowing only two out of the three parameters a, r and c. For 
example the Sierpinski Gasket has r — 3 and a — \/E, the T— fractal has r — 3 and a = y/E. 



IV. RANDOM WALKS AND HARMONIC OSCILLATIONS ON NON EXACTLY DECIMABLE 
FRACTALS: NT D AND THE CUTTING-DECIMATION RENORMALIZATION TRANSFORM 

In the previous section we have analyzed spectral dimension degeneration on exactly decimable fractals. Now we 
will consider the case of a class of fractal graphs (NTd) which are invariant under a more complex transformation 
T = D ■ C consisting of the product of a cutting transform C and a decimation D. We call T a cutting-decimation 
transform. It will be shown that Pa (e) and p{ui) behave differently under T and that NTjj are an explicit example of 
fractals with different diffusive and vibrational spectral dimensions. 

The fractal trees known as NTu || can be recursively defined as follows: an origin point O (Fig. 2) is connected to a 
point 1 by a link, of unitary length; from 1, the tree splits in k branches of length 2 (i.e. consisting of two consecutive 
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links); the ends of these branches split in k branches of length 4 and so on; each endpoint of a branch of length 2™ 
splits in k branches of length 2™ +1 . NTd can be naturally embedded in a suitable Euclidean space in such a way that 
their fractal dimension dp coincides with their connectivity dimension dc = 1 + lnfc/ln2 [ ppj . 

As one can easily verify, NTd are not exactly decimable since after a simple decimation starting from the origin O, 
one obtains k copies of the original structure joined together in a point instead of the same NTd- The NTd invariance 
under a T transform can be described in the following way (Fig. 3). Suppose to cut the log of the tree in point 1 
and to separate the k branches (cutting transform). Now, each branch can be obtained from the initial NTd by a 
dilatation with a factor 2. Eliminating all branches but one and decimating it (decimation transform), one obtains 
the original NTd- The T transform can be applied to solve the random walks problem. The cutting transform gives a 
relation between random walks on the whole tree and random walks on one of its branches; more precisely one relates 
^Oo e ( e )j t ne generating function of the probability of returning to point O after a random walk on the NTd tree, 
and p^ anch (e) , the generating function of the probability of returning to the starting point 1 after a random walk on 
one of the branches. This relation is given by jTo| |: 

Pootf = ~% — r 1 ^- ( 30 ) 

2ePft anch (e) + k 

The decimation transform on the surviving branch is obtained as for exactly decimable fractals by the transformation 
laws: 

e — * a 2 e = 4e 



por 



'" " nch (e) P'57" C V) = -P^ anch (e) = ~P$ anch (e) (31) 
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Now since the decimation procedure transforms a branch into the initial tree we have P' n (e') = Poo e ( e ') an< ^ 



relation (30) becomes 



P oo(e) = 2P £T:\ +k 02) 



4eP£g e (4e) + k 

Choosing a suitable expression for Pqq s {() |l0| we obtain from (p2|): 

,W + j|f (33) 

As for harmonic oscillations, the cutting transform gives a relation between p tree (u>) and p branch (a;) ; this can be 
obtained using the following properties pT] : 1) cutting or adding a finite number of points or bonds from an infinite 
network does not affect its vibrational spectrum; 2) the spectral density p(u>) (normalized to 1) of k copies of a given 
structure all attached in a point coincides with the corresponding one for a single structure. From these properties it 
immediately follows that the cutting transform does not affect the density of modes p(oj) so that: 

p tree (uj) =p branch (uj) (34) 

Applying the decimation transform we get a = 2 from ([3l]) while the decimation ratio is r = 2 so that ( pp[ ) becomes: 

p branch (to) -> p' branch {u') - p hranch (w) (35) 

Since p^ ranch {uo') = p tree (uj') relations (||) and (f|) give us p tree (uj) = p tree (uj') so that d v = 1. 
Since do = 1 + log kj log 2 and dy = 1, ATTd trees represent an explicit case of dynamical dimension splitting. All 
our results can be generalized |l2j to NTd lattices with a branches growth factor r (or decimation ratio) different 
from 2. In all these cases we will find again dynamical dimension splitting since dy — 1 and: 

d D = dy + ^ (36) 

logr 
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V. CONCLUSIONS 



In this paper we have shown how to built an exact renormalization group technique in real space using a more general 
reduction procedure than the usual decimation transform. This new technique allows to point out the phenomenon 
of dynamical dimension splitting which, as we have explicitly shown, is always absent in the case of exactly decimablc 
fractals, where do = dy < 2. Therefore these structures have very peculiar statistical mechanical properties and their 
study does not allow to distinguish between average and local quantities. 

Many examples of spectral dimension splitting can be found even on non fractal graphs, such as comb lattices 
and other branched structures. In general the situation with djj ^ dy is the most common so that it becomes 
particularly important to study the consequences of dynamical dimension splitting in all the phenomena which are 
influenced by the value of the spectral dimension such as diffusion, vibrational dynamics and phase transitions. 
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Figure captions 



Fig.l Exactly decimable fractals: a) Sierpinski Gasket, b) T— fractal c) Branched Koch curves 
Fig.2 NT D with k = 3 

Fig. 3 Cutting-Decimation procedure: 

a) Cutting of the log of the NT D 

b) Separation of the k branches 

c) Decimation of the points labelled by X 

d) Recovering of the original NTd 
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